Detection of entanglement with few local measurements 
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We introduce a general method for the experimental detec- 
tion of entanglement by performing only few local measure- 
ments, assuming some prior knowledge of the density matrix. 
The idea is based on the minimal decomposition of witness 
operators into a pseudo-mixture of local operators. We dis- 
cuss an experimentally relevant case of two qubits, and show 
an example how bound entanglement can be detected with 
few local measurements. 
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A central aim in the physics of quantum informa- 
tion is to create and detect entanglement - the resource 
that allows to realize various quantum protocols. Re- 
cently, much progress has been achieved experimentally 
in creating entangled states [1] . In every real experiment 
noise and imperfections are present so that the generated 
states, although intended to be entangled, may in fact be 
separable. Therefore, it is important to find efficient ex- 
perimental methods to test whether a given imperfect 
state p is indeed entangled. 

Obviously, the ultimate goal of entanglement detection 
is to characterize entanglement quantitatively, and iden- 
tify regions in the parameter space which allow to maxi- 
mize entanglement for a particular quantum information 
processing task. The first step towards this ambitious 
goal is to detect whether a given state is entangled or 
not. 

The question of direct detection of quantum entangle- 
ment has been recently addressed in Refs. [2-4]. In [3,4] 
the authors study the case of mixed states and find effi- 
cient ways to estimate the entanglement of an unknown 
state. Their method is based on structural approxima- 
tions of some linear maps followed by a spectrum esti- 
mation. Although experimentally viable the method is 
not very easy to implement and it requires further modi- 
fications in order to be performed by local measurements 
[5]. Here, we approach the same problem from a differ- 
ent perspective. We use special observables, the so-called 
witness operators [6,7] and their optimal decomposition 
into a sum of local projectors. Note that in this way 
we answer an open question posed recently in [8], where 
non-local measurements of entanglement witnesses were 
studied. 

The construction of a witness for a given arbitrary 
state is, in general, a formidable task. It can, however, 
be accomplished in typical experimental situations where 
one has some a priori information about the density ma- 



trix. This is always the case when the experiment is 
aimed at producing a certain state, rather than checking 
properties of an a priori unknown state. We discuss two 
experimentally relevant situations in this paper, namely 
the generation of a definite pure entangled state of two 
parties, and the generation of a specific bound entangled 
edge state. In both cases our method can be applied in 
arbitrary dimensions. 

Having constructed a witness, its measurement can be 
performed locally, since every observable can be decom- 
posed in terms of a product basis in the operator space. 
Here we propose two ways of optimizing such local mea- 
surements. The first one consists in looking for the opti- 
mal number of local projectors (ONP). The second one 
consists in searching for the optimal number of settings 
of detecting devices (ONS). By a setting of the device 
of a single observer we understand here the choice of 
the local orthonormal basis in the corresponding Hilbcrt 
space. The device measures then simultaneously pro- 
jections onto the vectors belonging to the basis; the set 
of these projectors forms a complete set of commuting 
observables [9]. A setting of the devices for a pair of 
observers corresponds then to a correlated choice of in- 
dividual settings. 

Both optimization methods are formulated as the prob- 
lem of decomposing a given operator into a sum of projec- 
tors on product states with an optimal number of terms. 
No general solutions for this kind of problems are known 
so far. 

Before describing the details of our method, let us 
briefly discuss other methods of entanglement detection 
with local measurements. For systems of two qubits or 
of one qubit and one qutrit, a necessary and sufficient 
criterion for entanglement, namely the non-positivity of 
the partial transpose [10], is known. Thus, using tomo- 
graphy of the state p, which can be achieved with local 
measurements [11], one can fully determine p, calculate 
the partial transpose and check its positivity. However, 
for two qubits this approach requires 9 different settings 
of the measuring devices in order to determine 15 pa- 
rameters describing the state in general. In our example 
below, only 3 settings suffice to detect entanglement if 
we have certain knowledge about p and we optimize the 
local decomposition of the corresponding witness. Our 
approach has similar advantages with respect to the de- 
tection of entanglement visibility [12], which in principle 
requires a continuous family of devices' settings. 

Another way of detecting entanglement by local mea- 
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surements consists in a test of a generalized Bell inequal- 
ity [13,14], with which our approach has a formal similar- 
ity. Nevertheless there exist many entangled states which 
do not violate any known Bell inequality [15]. However, 
for any given state one can always find a witness operator 
and its local decomposition, such that the entanglement 
of this state can be detected locally. For the situations, 
in which a previous knowledge of p allows to construct 
a suitable witness, our method is more powerful than a 
test of a Bell inequality. 

We introduce and illustrate our method in the sce- 
nario of creating and measuring the entanglement of 
two qubits, and then discuss shortly the generalisation 
to higher-dimensional states, including bound entangled 
states. 

Let us consider an experiment that produces the fol- 
lowing convex combination of a desired pure entangled 
state \ip}(ip\ and a mixed state a representing some noise, 



g = p\ i/>}(ip\ + (1 -p)a , 0<p<l 



(1) 



where | tp) can be written in the Schmidt decomposition 
as | ip) = a\ 01) + b\ 10) with a, b > and a 2 + b 2 = 1. 
The noise a is assumed to be within a ball of radius 
d around the totally mixed state, i.e. ||cr — 1/4-11 < d. 
Here ||A|| := y/Tr(A^A) is the Hilbert-Schmidt norm for 
operators A on the Hilbert space. One neither knows the 
probability p, nor the exact shape of a. The task is to 
determine whether g is entangled or not. 

Let us briefly summarise the well-known concept of en- 
tanglement witnesses [6,7]: a density matrix g is entan- 
gled iff there exists a Hermitian operator W such that 
Tr(Wg) < 0, but for all separable states Tr(Wg sep ) > 
holds. In this sense W "detects" the entanglement of 
g. Note that W has at least one negative eigenvalue. 
Methods to construct entanglement witnesses have been 
presented in Refs. [6,7]. For states with a non-positive 
partial transpose there is a simple and straightforward 
construction of W: let | e_) be the eigenvector of g TA that 
corresponds to its minimal (negative) eigenvalue, namely 
g TA \ e_) = X m in \ e-), with A mm < 0. Here Ta refers to 
partial transposition with respect to the first subsystem. 
Thus, W = (| e_)(e_ \) Ta detects the entanglement of g, 
as Tr((|e_)(e_ \) TA g) = Tr(|e_)( e _ \ Q Ta ) = X mm < 0. 
This witness is tangent to the set of separable states and 
already optimal [7], i.e. there is no witness that detects 
other states in addition to the ones detected by W. 

For g given in Eq. (1) and the case d = one finds 
Xmin = (1 -p)/4 -pab, and | e_) = ^-(| 00) - 1 11)), i.e. 
the witness is given by 



W = 



1 \ 
-10 
0-100 
1 / 



(2) 



Note that this witness neither depends on p, nor on a. 
Hence as long as the experimental apparatus produces 



any superposition of | 01) and 1 10) plus white noise, W 
will be a suitable operator. For d ^ it still provides the 
possibility of entanglement detection. 

Let us point out that this witness is also suitable for 
other physical scenarios. For example, consider the case 
where the noise mechanisms in the experimental setup 
are characterised by memory effects, or the case where 
the entangled state \tp), (defined via Eq. (1)), is gen- 
erated perfectly and then sent through a transmission 
channel with correlated noise. If the noise mechanisms 
acting on the state can be described as a depolarising 
channel with some correlations of strength /j, [16], then 
the resulting state will be of the form: 



g = {1® 1 + ??(a 2 - b 2 )[a z <g> 1 - 1 1 

+ [m+(1- v)v 2 ][-(Tz ® <?z 

+ 2ab(a x ®<j x +<j v ® a y )]} /4 . 



(3) 



Here rj and \i describe the depolarisation and the degree 
of memory introduced by the noise process. This fam- 
ily of states is now characterized by three independent 
parameters: a, r] and p. Remarkably, W turns out to 
detect the entanglement of this whole family of states. 
This is proven by calculating the range of the family's 
three parameters where g TA has a negative eigenvalue, 
and showing that for this range Tr(Wg) < holds. 

We look now for a decomposition of the witness into a 
sum of projectors onto product vectors, i.e. 

W = ^2c l \a i ,b i )(a i ,b l \ = ^ a\ a^fa | <g> | 6i><6» | , (4) 



where the coefficients Cj are real and fulfill Cj = 1 . 
Note that at least one coefficient has to be negative - 
this characterises a so-called pseudo-mixture. Any bipar- 
tite Hermitian operator can be decomposed in projec- 
tors onto product states, like in Eq. (4), in many dif- 
ferent ways. However, we are interested in finding the 
optimal decompositions in the two ways described above. 
Optimal pseudo-mixtures in the sense of minimizing the 
number of non- vanishing q represent an ONP. They have 
been studied for two qubits in [17] where it was shown 
that any general vector \<f>) — a\ 00) + (3\ 11) with a, (3 
real and different from zero, and a 2 + 1 — 1, can be 
decomposed minimally with 5 terms: 



i=l 

-a/3(|01)(01 | 4 
where we have used the definitions 



io><io |) 



(5) 



|/i) =e-^ CO s0|O)+e 4 * sin0| 1) = | f 2 *) 

|/ 3 >=cos0|O>+sin0|l> = |/i> + l/2> 
cos6 = y/a/{a + (3), sine = ^J(3/{a + f3). (6) 

Here * denotes complex conjugation. The case a = — 
—(3 corresponds exactly to the decomposition of W in 
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Eq. (2) that we are looking for. Such a decompostion 
into 5 terms requires four different settings of the mea- 
suring device (in the sense described above and in the 
footnote [9]). It is therefore optimal in the number of 
projectors (ONP), but not optimal with respect to the 
number of correlated devices' settings (ONS). Indeed, the 
witness W considered in Eq. (2) can be optimally imple- 
mented by using only three settings. This can be shown 
as follows: defining the eigenstates of the Pauli matri- 
ces as | z+) = | 0), | z~) = | 1), | x±) = ^(| 0) ± 1 1)) and 

|» ± ) = ^(|0)±i|l)) ) we find 

(| 0) (0 |) T -4 = a 2 \ z + z+)(z+z + | + f3 2 \ z~z-)(z-z- | 
+a/3(| x + x + ) (x + x + | + | x~x~) (x~ x~ \ 

-\y + y~)(y + y~ \ - \y~y + )(y~y + \) ■ (?) 

Note that this decomposition contains 6 terms, but only 
three correlated devices' settings. Alice and Bob have to 
classically correlate their measurements in their respec- 
tive x, y- and z-directions as indicated in (7), and add the 
resulting expectation values with the according positive 
or negative weight in order to determine Tr(Wg). Let us 
not in passing that in general an entanglement witness 
does not provide an entanglement measure. However, for 
fixed ab and d = the probability p in Eq. (1) can be 
found from the measurement outcome via the relation 
p= {l-4Tr(Wg))/(l + 4ab). 

The pseudo-mixture (7) provides an ONS, since it is 
impossible to decompose (| 4>){<j) \ ) Ta with less than three 
devices' settings. Let us assume the contrary, i.e. 



{\m\) TA = Y, u ij \u?){u?\®\ u f){uf\ + 

i2^\v?)(v?\®\vf){vf I, (8) 



where (i 



l ) = S a = (t 



l ) and the same holds for 



B. We expand (\4>) (4>\) Ta 
we denote er = 1, and 



=0 ^i] 



o~i (g> <jj , where 
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2 u 

o h 



(9) 
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Note that the 3x3 submatrix in the right bottom 
corner is of rank three. Now we write any projec- 
tor in the rhs of (8) with Bloch vectors: the pro- 
jector (wfXwfl = Y^\=o s t a i IS represented by the 
vector s A = |(1, sf-, s% , s£ ), and by sj_ = 

|(1, — Si, — sf, — s^). Expanding the first sum on the 
rhs of (8) in the (<7, <g> aj ) basis leads to a 3 x 3 subma- 
trix in the right bottom corner which is proportional to 
((7ii-C/i 2 -C/2i + C/22)(s^,s^,^) T ( S f, S f, S f),and thus 
of rank one. The corresponding matrix from the second 



sum on the rhs of (8) is also of rank one, and we arrive at 
a contradiction: no matrix of rank three can be written 
as a sum of two matrices of rank one. □ 
We now emphasize the power of witness operators as a 
tool for the detection of entanglement by discussing the 
noise in Eq. (1) in some more detail. When d ^ 0, the 
state g lies within a ball B Pt d with radius (1 — p)d. If 
p is such that this ball is either included in the set of 
separable states, or in the set of entangled states, then 
the given W is optimal and the sign of Tr(Wg) provides 
a signature of entanglement versus separability. If, how- 
ever, the ball lies across the boundary between those two 
sets, errors may occur. Different questions can then be 
addressed: first, one may want to be sure that a given 
state is separable. For the case a = b = l/s/2 (which 
we assume here and in the following) we can estimate a 
lower bound t such that if Tr(Wg) > t then g(p, d) is 
necessarily separable. This bound, which depends on d, 
is given by: 



r(d) = \-<?-^j{^-<P){\-d?). (10) 

For any r' with < r' < r, there exists an entangled 
state g(p,d) with Tr(Wg) = r' . To derive Eq. (10), one 
uses the fact that there is a ball B of separable states of 
maximal radius 1 / \/l2 around 1/4. Since we do not know 
p we can only say that g e U p e[o 1] ^p,d> which is a kind 
of a convex cone originating in | ip) (i/j | and terminating in 
the ball B o d of radius d around the totally mixed state. 
The bigger Tr(Wg), the closer to the ball Bo,d C B is g. 
Obviously, if Tr{Wg) is big enough we have g £ B. In 
this manner one can determine the value of r. □ 

Second, one may be interested in minimizing the prob- 
ability of making an error, either by mistaking a separa- 
ble state for an entangled one, or vice versa. If we assign 
Tr(Wg) > <=> g separable, it turns out that - depend- 
ing on the value of d - in order to minimize this error, it is 
more favorable to use W e := W—el. This operator is not 
a "witness" in the original sense, because it yields nega- 
tive expectation values for some separable states. Note 
that it requires the same measurement settings as W. 
To estimate the error inherent in this detection scheme, 
we have used the method from [18] to randomly gener- 
ate a sample of 50000 density matrices of the form (1), 
and then checked their separability using both the partial 
transposition criterion and applying W £ . The percentage 
of errors when using W e is plotted in Fig. 1. For large d 
the operators W e are in fact less erroneous in detecting 
entanglement. Further numerical analysis suggests that 
the optimal e increases quadratically with d. 

Generalizing the above results to higher dimensions for 
states with non-positive partial transpose (NPT states) 
is possible, although not straightforward. In an N x M 
dimensional Hilbert space with N < M, we first have to 
identify the vector \<p) that corresponds to the minimal 
(negative) eigenvalue of g TA . Without loosing generality 
we may assume that it has maximal Schmidt rank, and 
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FIG. 1. Error (maximised over all p) as a function of e, for 
three different values of d. 

is given by \4>) = X^o, 1 a i\^)- Obviously, the ONP cor- 
responding to i\(t>){4>\) must contain at least TV 2 terms 
since the rank of (\<f>) (<P\) Ta is TV 2 . For TV = 2, i.e. for 
2x M dimensional systems the results obtained for M = 2 
are also valid for M > 2, since the maximal Schmidt rank 
in such spaces is 2. This implies that the ONP must con- 
tain 5 terms, whereas the ONS can be realized with 3 
settings. For TV x M systems with TV > 3, TV < M 
we can easily construct a pseudo-mixture with 27V 2 — TV 
terms using the same method as in the case of 2 x 2 
systems. This gives the upper bound for the number of 
terms in ONP, and corresponds to an upper bound for 
ONS of 27V - 1 (27V) for even (odd) TV. By generalizing 
the method used to demonstrate that two settings are 
not enough in 2 x 2 to the TV x M case one can prove 
that any ONS must contain at least TV + 1 settings. It 
is not clear, however, whether this bound can be reached 
in general. 

In higher dimensions there also exist entangled states 
with positive partial transpose, namely bound entangled 
states [19,20]. For this type of states no general oper- 
ational entanglement criterion is known, and thus even 
the full knowledge of the density matrix may not suf- 
fice to decide whether the density matrix is entangled or 
not. There exists, however, an important class of bound 
entangled states, the so-called "edge" states, for which 
optimal witness operators can be constructed explicitly. 
A state S is called an edge state iff it cannot be rep- 
resented as 5 = qS' + (1 — q)a, where a is a separable 
state, 5' is a state with positive partial transpose, and 
< q < 1. Edge states are, in a certain sense, the bound 
entangled analogues of pure entangled states. In the sit- 
uation where an experiment is aimed at the generation 
of an edge state, our method of local decomposition of a 
witness provides a genuine experimental test. 

Let us illustrate this with the example of unextendiblc 
product basis (UPB) states [21] in a 3 x 3 dimensional 
space. The states 



|^o) = _|0)(|0)-|l)), 



|V> 2 ) = - 7 =|2)(|1)-|2)), 



|Vi) = -^(|0)-|l))|2), |V> 3 ) = ^=(|1)-|2))|0), 



|V>4> = 3(|0> + |1> + |2»(|0> + |1> + |2» 



(11) 



form a UPB, i.e. they are orthogonal to each other and 
there exists no product vector which is orthogonal to all 
of them. The state p BE = j(l - Yh=o I V>i)(^ I) con- 
structed from this UPB is an entangled state with posi- 
tive partial transpose. The generic form of an entangle- 
ment witness for such a state is [6,71 



W= (P• 



inf<e,/|P 

I e./> 



TA )/2-et 

Q Ta 



(12) 



e,/)/2 , 



where P and Q denote the projectors onto the kernel of 
Pbe and the kernel of p T B A E , respectively. For the given 

UPB state we have P = Q Ta = YlUo I (V'i I ■ The main 
problem for the construction of W is to find e. An analyt- 
ical bound obtained by Terhal [6] gives e > 0.0013. Nu- 
merical analysis leads however to the much bigger value 
e ~ 0.0284. Once e is found, the decomposition of W is 
straightforward: the explicit form of the witness fixes five 
elementary measurements, and the identity can be de- 
composed into nine orthogonal projectors onto product 
vectors such that four of them coincide with four of the 
UPB states. This follows from the fact that the vectors 
\i>o): IV'i)) 1^2), |^3) can be extended to an orthonormal 
basis by defining 



|^ 4 > = _|0>(|0> + |1», 



1 



V2 

&> = |1>|1>. 



(|0> + |1»|2>, 



1^5) 
1^7) 



V2 
1 

71 



|2>(|1> + |2», 

(|1) + |2))|0), 
(13) 



Altogether we are left with a pseudo-mixture that con- 
tains 10 projectors. Denoting by P>\ = {| 0), 1 1), | 2)}, 
B2 = {(|0> - |1»/V2,|2},(|0} + |1})/V2}, B 3 = 
{(| 1) - I 2))/V2, 1 0), (| 1) + I 2)/V2)}, and B 4 = {(| 0) - 
I 1»/V2, (| 0) + I 1) + I 2»/V3, (| 0) + 1 1) - 2| 2»/2}, we 
easily see that measurement of W for this decomposition 
requires 6 correlated settings for Alice and Bob: BiB 2: 
B 2 B U BiB 3 , B 3 Bi, B A B A , and B X B X . This result im- 
plies that the ONS must be < 6. By subtracting in Eq. 
(12) some positive operator / of full rank instead of the 
identity 1, one can reduce the number of projectors in 
the decomposition of W to 9, and this gives an ONP, 
since the number of terms in any ONP must be larger or 
equal than the rank of the witness, which is equal to 9. 
The idea is to form J as a convex sum of projectors onto 
\ipi)i=o,...,4, an( i onto 4 other product vectors that are ob- 
viously not orthogonal to |V'i}i=o,...,4) but such that the 
set of the 9 vectors forms a basis. Here the bound for e 



4 



has to be adapted to the choice of /, such that positivity 
of W on all separable states is guaranteed, i.e. 



inf 

\ej) 



(e,f\(P + Q T -)/2\e,f) 
(e,f\I\e,f) 



(14) 



If we choose as additional vectors |V>i)i=4,...,7 the de- 
composition contains 9 projectors in only 5 settings. Nu- 
merical analysis leads to e' ~ 0.0311 for this choice of 
I. Note that when the bound entangled state is affected 
by white noise, namely p p = p ■ pbe + (1 — p)l/9, the 
witness given above is still suitable for the detection of 
entanglement, provided Tr(Wp p ) < 0. For the witness 
in Eq. (12) this is the case when p > (1 — 9e/5). Let 
us mention that for experimental purposes it is not nec- 
essary to decompose the identity in Eq. (12) since this 
term will only add a constant (— e) to the outcome of 
the measurement of the "prewitness" W = (P + Q Ta )/2. 
This implies that indeed one only needs to find the corre- 
sponding ONP and ONS for W. In this way one obtains 
an ONP with only 5 terms, and that the ONS has to be 
less or equal to 5. 

In summary, we have introduced optimal decomposi- 
tions of witness operators into local projectors for the de- 
tection of entanglement. This method can be used with 
present experimental techniques. It is a very powerful 
method to detect entanglement in the cases where one 
has a certain knowledge about the state that one wants 
to create, e.g. when the aim is to produce a specific pure 
entangled state, but this state is corrupted by noise. At 
the present stage of quantum information processing, sev- 
eral experiments strive at creating such entangled states, 
and thus it is important to show that the produced state 
is indeed entangled. The more knowledge about the state 
is given, the less knowledge about the underlying noise is 
necessary for unambiguous classification of the state. For 
the situation in which little knowledge about the created 
state exists it will be favourable and maybe necessary to 
utilize more than one witness operator. A detailed anal- 
ysis of the trade-off between the initial knowledge of the 
state and the witnesses needed is left for further research. 

After submission we became aware of a recent preprint 
by A. Pittenger and M. Rubin [22] where the ideas of 
this paper have been further developed. In particular 
it is shown there that if N is prime, a projector onto a 
maximally entangled state with full Schmidt rank can be 
measured with N + 1 local measurements, so our bound 
can be reached for this case. 
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